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ABSTRACT 

Some supernovae (SNe) show evidence for mass-loss events taking place prior to their explosions. 
Measuring their pre-outburst mass-loss rates provide essential information regarding the mechanisms 
that are responsible for these events. Here we present XMM-Newton and Swift X-ray observations 
taken after the latest, and presumably the final, outburst of SN2009ip. We use these observations as 
well as new near infra-red and visible light spectra, and published radio and visible light observations 
to put six independent order-of-magnitude constrains on the mass-loss rate of the SN progenitor 
prior to the explosion. Our methods utilize: the X-ray luminosity, the bound-free absorption, the 
Ha luminosity, the SN rise-time, free-free absorption, and the bolometric luminosity of the outburst 
detected prior to the explosion. Assuming spherical mass-loss with a wind density profile, we estimate 
that the effective mass-loss rate from the progenitor was between 10 -3 to lO _2 M0yr _1 , over a few 
years prior to the explosion, with a velocity of ~ 10 3 kms -1 . This mass- loss rate corresponds to a total 
circum stellar matter mass of ~ 0.04 M , within 6 x 10 15 cm of the SN. We note that the mass-loss 
rate estimate based on the Ha luminosity is higher by an order of magnitude. This can be explained 
if the narrow line Ha component is generated at radii larger than the shock radius, or if the CSM has 
an aspherical geometry. We discuss simple geometries which are consistent with our results. 
Subject headings: stars: mass-loss — supernovae: general — supernovae: individual (SN2009ip) 



1. INTRODUCTION 

Supernova (SN) observations, especially of Type Iln 
(e.g., Filippenko 1997), indicate that some massive stars 
lose considerable amounts of mass (<; 10 _4 M Q ) within a 
few months to years prior to their explosions (e.g., Dopita 
et al. 1984; Chugai et al. 1994, 2004; Ofek et al. 2007, 
2010, 2013b; Smith et al. 2007, 2008, 2009; Kiewe et 
al. 2012). Several theoretical mechanisms to eject large 
amounts of mass with super-Eddington luminosities have 
been suggested. Quataert & Shiode (2012) suggest that 
in some massive stars the super-Eddington fusion lumi- 
nosities, shortly prior to core collapse, can drive con- 
vective motions, that in turn excite gravity waves that 
propagate toward the stellar surface. The dissipation 
of these waves can unbind up to several solar masses 
of the stellar envelope. In Ofek et al. (2013b) we ar- 
gued that this mechanism can unbind a lower amount of 
mass (~ 10" 2 M o ). Arnett & Meakin (2011) suggested 
that shell oxygen burning in massive stars produces large 
fluctuations in the turbulent kinetic energy, that in turn 
may produce bursts. Chevalier (2012) suggested that the 
mass loss is driven by a common-envelope phase due to 
the inspiral of a neutron star into a giant companion core, 
unbinding the companion envelope and setting up accre- 
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tion onto the neutron star, that in turn collapses into a 
black hole and triggers a SN explosion. Soker & Kashi 
(2013) suggested that the SN 2009ip explosion was due to 
the merger of two stars, while some of the pre-explosion 
outbursts occurred near periastron passages of the binary 
system. Another possible mechanism is the pulsational 
pair instability which in very massive stars can generate 
several explosions, expelling ^ 1 M each, followed by 
the collapse of the stellar core (Rakavy, Shaviv & Zina- 
mon 1967; Woosley, Blinnikov & Heger 2007; Waldman 
2008). 

Measuring the mass-loss rates from massive stars prior 
to their explosion can be used as a tool to study the latest 
stages of stellar evolution, and to discriminate between 
the different models suggested to generate large mass- 
loss events. Objects in which super-Eddington outbursts 
were directly observed prior to the SN explosion provide 
a way to constrain the time at which mass-loss was tak- 
ing place, and relate the optical luminosities with mass- 
loss rates and kinetic energy estimates. To date there 
are only three SNe in which precursor outbursts were 
seen prior to the SN explosion. These are the Type Ibn 
SN2006jc (e.g., Foley et al. 2007; Pastorello et al. 2008), 
the Type Iln SN2009ip (e.g., Mauerhan et al. 2012; Pa- 
storello et al. 2012; Prieto et al. 2012) and the Type Iln 
SN2010mc/PTF10tel (Ofek et al. 2013b). 

Here we present XMM-Newton and Swift X-ray obser- 
vations of SN 2009ip. We use these observations as well 
as published and new visible light and radio observations, 
to set an order of magnitude estimate on the mass-loss 
prior to the SN explosion. 

SN 2009ip was a Luminous Blue Variable (LBV) orig- 
inally detected in outburst on 2009 August 26.11 by the 
CHASE survey at a projected distance of 4.3 kpc from 
NGC7259 (Maza et al. 2009). Three additional out- 
bursts were subsequently discovered with the Catalina 
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Real-Time Transient Survey on 2010 July 15, on 2010 
September 29 (Drake et al. 2010), and then again on 
2012 July 24 (Drake et al. 2012). Based on its multiple 
outbursts, Smith et al. (2010) and Foley et al. (2011) 
argued that it is a supernova impostor (see recent re- 
views in Kochanek, Szczygiel, & Stanek 2012; Smith et 
al. 2011; van Dyk & Matheson 2012). In September 
2012, Smith & Mauerhan (2012), and later Mauerhan et 
al. (2012), reported the detection of broad P-Cygni lines 
with velocities of up to 13,000 km s _1 , suggesting that 
the star had finally exploded as a Type Iln SN. Previ- 
ous cases in which a likely LBV progenitor has exploded 
as a supernova include SN 2005gl (Gal- Yam et al. 2007; 
Gal- Yam & Leonard 2009), and SN 1961V (Kochanek, 
Szczygiel, & Stanek 2011; Smith et al. 2011). Prieto et 
al. (2012) reported that around 2012 September 24 the 
object's /-band light curve started to rise rapidly at a 
rate of 2.3magday~ 1 . Shortly afterwards, on early Oc- 
tober, X-ray emission was detected from SN 2009ip with 
the Swift/X-Ray Telescope (XRT; Margutti & Soderberg 
2012b). 

Throughout the paper we assume that the source is 
located at a distance of 20.4 Mpc. In Sj2]we present our 
observations of SN2009ip, while in $3] we review various 
methods for estimating the mass content of the SN cir- 
cumstellar matter (CSM). Finally, in <j4] we apply these 
methods to SN 2009ip and discuss our findings. 

2. OBSERVATIONS 

We observed SN2009ip with VMM-Newton on 2012 
November 1 in prime full window imaging mode for an 
effective exposure time of 8 ks. Using data collected with 
the EPIC-pn detector (Tian et al. 2007), we accumu- 
lated the source spectrum from a circular region of 30" 
centered on the optical position of SN2009ip. We se- 
lected a circular background region from a source free 
area on the same chip (i.e., CCD 7) with the same aper- 
ture size. The source is detected at a significance of 
about 3 a with a background subtracted count rate of 
(4.6±1.5)xl0~ 3 counts s _1 , yielding a total of 37 net 
source counts in the 0.5—10 keV. We generated the de- 
tector and ancillary response files using the latest cali- 
bration data. 

The Swift-XRT (Gehrels et al. 2004) observed 
SN 2009ip on an almost daily basis since 2012 September 
4 (triggered by Roming/Maragutti). Some of these X-ray 
observations have been already reported in e.g., Margutti 
et al. (2012), and Campana (2012). For each Swift- 
XRT image of the SN, we extracted the number of X-ray 
counts in the 0.2-10 keV band within an aperture of 9" 
radius centered on the SN position. We note that this 
aperture contains ks 50% of the source flux (Moretti et al. 
2004) . The background count rates were estimated in an- 
nuls around the SN location, with an inner (outer) radius 
of 50" (100"). The log of Swift-XRT observations, along 
with the source and background X-ray counts in the in- 
dividual observations are listed in Table [T] SN 2009ip 
is only marginally detected in individual images, but it 
is clearly visible in the coadded data. Figure [1] shows a 
binned light curve based on the Swift-XRT observations. 
The binned measurements are listed in Table [5] 

For our XRT spectral analysis we selected all the XRT 
observations between MJD 56174 and 56228, taken in 
photon counting mode and with an integration longer 
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Fig. 1.— X-ray light curve of SN 2009ip based on Swift-XRT 
(circles) and XMM-Newton (square) observations. The triangles 
mark XRT 2-cr upper limit. The horizontal error bars represent 
the range of observations in each bin. The gray dashed line in- 
dicates the mean XRT count rates level of the observations taken 
between 2012 September 29 and 2012 Nov 28. It corresponds to 
a luminosity of « 1.1 X 10 39 ergs -1 . We note that the left-hand 
axis count rate corresponds only to the Swift-XRT observations. 
The right-hand axis shows the unabsorbed luminosity assuming, a 
Galactic hydrogen column density of Nu = 1.2 X 10 20 cm -2 , and 
an X-ray spectrum of the form n(E) <x E~ 1S , where n(E) is the 
photon numbers per unit energy. The black solid line shows the 
expected, order of magnitude, evolution of the X-ray luminosity 
assuming optically thin wind-profile CSM with mass-loss ra te o f 
7 X 1Q-*M Q yr" 1 and v w = 500 kms -1 (based on Eqs.l8land[j"2l). 



TABLE 1 

Swift-XRT OBSERVATIONS 



MJD 


Exposure time 


Source 


Background 


day 


ks 


ct 


ct 


55084.44 


9.86 





19 


56174.86 


1.96 





1 


56176.60 


1.78 





4 


56183.40 


1.65 





4 


56190.75 


0.39 









Note. — MJD is the modified Julian day. Source is the number 
of counts in a 9"-radius aperture of the source position and in the 
0.2-10 keV band. Background, is the number of counts, in the 0.2- 
lOkeV band, in an annuls of inner (outer) radius of 50" (100") 
around the source. The ratio between the background annulus 
area and the aperture area is 92.59. This table is published in its 
entirety in the electronic edition of ApJ. A portion of the full table 
is shown here for guidance regarding its form and content. 

than 500 s. This resulted in a total effective exposure 
time of 149 ks. We extracted a stacked source spec- 
trum from a circular region centered at the SN location 
with a radius of 20". The stacked background spec- 
trum was accumulated from a 20" circular source-free 
region for all observations. We grouped the source spec- 
trum with a minimum of 10 counts in each energy bin. 
The background subtracted count rate of the source is 
(4.1 ± 0.7) x 10~ 4 counts s , corresponding to 61 net 
source counts. 

We used XSPECfl V12.7.1 (Schafer 1991) to simultane- 
ously fit the VMM and Swift-XRT spectra. In all cases 
we set the Galactic extinction to Nn = 1.2 x 10 20 cm~ 2 

7 http:/ /heasarc. gsfc.nasa.gov/xanadu/xspec/ 
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Fig. 2. — Left: Ay 2 contours for fitting the XMM and Swift-XKT spectra with a Mekal model in the column density vs. temperature 
(kT) space (see text for details). The contours represent the 1, 2 and 3-<r errors, while the plus sign represents the best fit model, middle: 
The same as the left panel, but for a power-law model. The X-axis represents the power law index, T. Right: The same as the left panel, 
but for a black-body model. The X-axis represents the black-body temperature in keV. 



TABLE 2 
Swift-XKT BINNED DATA 



(MJD) 


Ranj 


;e 


CR 


ul 2ct 


Exp. 


day 


day 


day 


cnt/ks 


cnt/ks 


ks 


55084.4 


-0.0 


0.0 




0.61 


9.86 


56178.0 


-3.2 


12.7 




1.03 


5.79 


56195.4 


-2.7 


2.0 




0.55 


17.40 


56201.6 


-3.9 


2.9 


<52+ - 20 




46.38 


56207.3 


-1.9 


2.9 


n qv+0.29 
u ' al -0.23 




37.17 


56214.5 


-2.3 


4.3 


35+ ' 28 
u "->°_0.17 




22.66 


56222.8 


-2.5 


5.9 


50+ ' 30 
U ' OU -0.20 




24.08 


56233.9 


-3.5 


10.5 


u.oo_ . 20 




30.80 


56253.3 


-6.4 


5.2 


70+ ' 48 

u - ,u -0.30 




14.20 


56279.7 


-19.4 


14.7 




0.18 


52.57 



Note. — Binned Swift-XKT light curve of SN2009ip. (MJD) is 
the weighted mean modified Julian day of all the observations in a 
given bin, where the observations are weighted by their exposure 
times. Range is the time range around (MJD) in which the light 
curve (Table [TJ was binned. CR is the counts rate along with the 
lower and upper 1-<t errors. We note that the source count rates 
are corrected for extraction aperture losses. UL2o- is the 2-cr upper 
limit on the source count rate, which is given if the total source 
counts within the binned exposure is < 1. Exp is the exposure 
time. 



TABLE 3 

X-RAY SPECTRAL PARAMETERS 



Model 




parameter 


cm * 


x 2 /d.o.f. 


Mekal 
Power law 
Black body 


kT 

r = 

kT 


= 4.74+2 8 3 keV 
i 7Q+0.60 
i -' M -0.50 

= 0.72 ±0.10 keV 


(2.8± 2 '*) x 10 21 
(3.2± 2 ;J) x 10 21 


14.66/13 
14.98/13 
18.19/13 



Note. — r is defined as the power-law index in a spectrum of 
the form, n(E) oc E~ r , where n(E) is the number of photons per 
unit energy, d.o.f. is the number of degrees of freedom. 



(Dickey & Lockman 1990), and wc fitted four parame- 
ters: normalization of the XMM spectrum, normaliza- 
tion of the Swift-XRT spectrum, a parameter describing 
the spectrum (i.e., temperature or power-law index), and 
the hydrogen column density at the redshift of the super- 
nova (z = 0.00594) assuming solar metallicity. The best 
fit parameters are listed in Table |3] Figure [5] 

shows the A% 2 contours of these fits in the Nn vs. tem- 
perature (kT) or power-law index space. From these fits 



we can set a 3-ct upper limit on Ah in the SN CSM of 
2xl0 22 cm- 2 . 

On 2012 September 22, two days before the fast rise in 
the light curve, we obtained a near-infrared (NIR) spec- 
trum of SN 2009ip with the Folded-port InfraRed Echel- 
lette spectrograph (FIRE; Simcoe et al. 2008, 2010) on 
the 6.5-m Magellan Baade Telescope. We used the low- 
dispersion, high-throughput prism mode, and completed 
an ABBA dither sequence. The data span 0.8-2.5 /im at 
a resolution ranging from 300-500. Immediately after- 
wards, we obtained a spectrum of an A0V standard star 
for the purposes of flux calibration and removal of telluric 
absorption features, as described by Vacca, Cushing, & 
Rayner (2003). Data were reduced using the FIREHOSE 
pipeline developed by R. Simcoe, J. Bochanski, and M. 
Matejek. Smith et al. (2013) present detailed analysis of 
the NIR spectrum. 

On 2012 December 4, we obtained a visible light spec- 
trum of SN 2009ip using the Dual Imaging Spectrograph 
(DIS) mounted on the ARC-3.5m telescope at Apache 
Point Observatory to get a spectrum, with 600 s integra- 
tion, of SN2009ip in the wavelength range from 3500A to 
9000 A and resolution of about 400. The visible-light 
spectrum was flux calibrated using the standard star 
BD+28°4211 Parts of the IR and visible light spectra, 
centered on the Paschena and Ha lines, respectively, 
are shown in Figure [31 The full spectra are available 
from the WISeREP archive (Yaron & Gal- Yam 2012). 
We fitted a two Gaussians model to the Paschena and 
Ha lines. We find that in the Paschena, the narrow 
(wide) component width corresponds to velocity of « 200 
(» 2100) kms , In the Ha line, the narrow line com- 
ponent width corresponds to velocity of ~ 300 kms , 
while the difference between the emission line center and 
the bottom of the P-Cygni absorption feature is about 
8000 km s _1 . The flux of the Ha narrow component is 
about 3 x 10~ 14 ergcm -2 s _1 . Since we do not have ac- 
cess to photometric measurements of the SN obtained 
around the same time in which we got the visible light 
spectrum, we estimate that the line flux measurement is 
good to about 30%. For future calibration, we note that, 
based on the current calibration, the SDSS AB synthetic 
magnitudes of the visible light spectrum are 17.17, 16.31 
and 16.11 in the g, r and i-bands, respectively. 

3. MASS LOSS ESTIMATORS 
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Fig. 3.— Left: Near IR spectrum of SN2009ip. The X-axis 
is the wavelength at the observer's frame. The spectrum is cen- 
tered on the Paschena line. The vertical lines on the scale cor- 
responds (from left to right) to velocities of — 10 4 , — 10 3 , 0, 10 3 
and 10 4 kms — 1 relative to the line center. The dashed-gray line 
shows the best-fit two-Gaussian model (we note that the Gaussian 
were convolved with the approximate instrumental broadening). 
The gaps in the plot are due to the removal of bad/noisy pixels. 
Right: Same as the left spectrum, but for the visible-light Ha line. 



Here we review several methods that we use to con- 
strain the mass-loss rate from SN progenitors. In i|4]we 
implement these methods for the case of SN 2009ip. 

We use several observables, including the bound-free 
absorption limit derived from the X-ray spectrum, the X- 
ray luminosity, upper limit on the diffusion time scale as 
derived from the SN rise time, the Ha luminosity, the non 
detection in radio bands, and the bolometric luminosity 
of the precursor observed prior to the SN explosion. 

Our model assumes that the interaction of the SN blast 
wave with the CSM produces X-ray and radio emission 
at the radius of the shock. The visible light photons may 
be produced below, or at, the shock radius. The Ha 
emission can be emitted either at the shock region, if it 
is due to collisional processes, or above the shock radius 
if it originates from optically thin CSM ionized by the SN 
radiation field. All the photon diffusion and attenuation 
is taking place above the shock. 

Throughout the paper we assume that the CSM around 
the progenitor has a spherical wind-density profile of the 
form p = Kr~ 2 , where r is the distance from the progeni- 
tor, K = M /(47ru w ) is the mass-loading parameter, M is 
the mass- loss rate, and u w is the wind/outburst velocity. 
Given the outbursts observed in SN2009ip prior to its 
final explosion, it is likely that the CSM was not ejected 
as a continuous wind with a uniform velocity. However, 
the mass-loss rate estimators we use below, are not very 
sensitive to this assumption. The reason for this is that, 
for a reasonable density distribution, the emission, or 
attenuation, are calculated by integrals which are dom- 
inated by the value at the shock radius. Here, the only 
exception is the mass-loss estimator based on the Ha line 
luminosity (see i i3.2[) . Therefore, we argue that the use 
of the continuous wind density profile provides an or- 
der of magnitude estimate for the mass-loss rate. In the 
following sections we discuss our specific mass-loss rate 
estimators and their caveats. 



3.1. Bound free absorption 

The particle density profile, in continuous wind, is 
given by (e.g., Chevalier 1982) 



n 



1 



M 



(li p ) 4irm p v w r 
(Mp) 



-1 -2 -3 
'w,500 r 15 cm ' 



(1) 



where Mrj.oi is the mass-loss rate in units of 
10~ 2 Mq yr _1 , w w ,500 is the wind/outburst speed in units 
of 500 km s -1 , ris is the radius in units of 10 15 cm, m p is 
the proton mass and (/i p ) is the mean number of nucleons 
per particle (mean molecular weight). For our order-of- 
magnitude calculation, we adopt (/i p ) = 0.6. In a wind 
profile, the column density between radius r and infinity 



N = 



(2) 



Assuming that the gas in the pre-shocked wind is neu- 
tral and has solar abundance, the bound-free optical 
depth in the 0.03-10 keV region is roughly given by (e.g., 
Behar et al. 2011)0 



-na = Na(E) 



'SQM 0M v-] 50Q r^El 



-2.5 



(3) 



where <j(E) is the bound- free cross section as a function 
of energy E, and E\ is the energy in keV. This approx- 
imation is valid when the material is neutral. However, 
since above ~ 0.5 keV, metals with a high ionization 
potential dominate the absorption, this formula is still 
valid, to an order of a magnitude, above 0.5 keV when 
some of the inner electrons of the metals are bound (i.e., 
even if all the hydrogen is ionized). Chevalier & Irwin 
(2012) estimated that the metals will be completely ion- 
ized only above shock velocities of about 10 4 kms _1 . 

3.2. Ha luminosity 

Assuming the SN radiation field can ionize all the hy- 
drogen in the CSM, the mass of the hydrogen generating 
the Ha line is 

Mh«^. (4) 

Here h is the Planck constant, Lh q is the Balmer Ha line 
luminosity, i/jj is the line frequency (4.57 x 10 Hz for 
Ha), and a| s is the case-B effective recombination coef- 
ficient at 10,000 K (w 8.7 x 10~ 14 cm 3 s" 1 for Ha; Oster- 
brock & Ferland 2006). An important caveat is that this 
estimate assumes that the line is generated by recom- 
bination. Any Ha radiation generated in the shocked 
region due to collisional processes, is not included here. 
In order to avoid this problem we take as the line lumi- 
nosity only the luminosity of the narrow line component, 
which we assume is due to wind above the shock region. 
In a wind profile, the integrated mass from radius r to 
ri is 



M: 



l-ri 

- / 4:nr 2 Kr~ 2 dr = AiiK{ri — r) ~ AirKfir 

J r 



8 This approximation deviates by a factor of two from a more 
accurate calculation (e.g., Morrison & McCammon 1983). 
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(5) 



where j3 = (n— r)/r. We note that f3 cannot be arbitrar- 
ily large (otherwise the mass in the CSM will diverge), 
and here we will assume it is of order unity. 

By substituting Equation [T] into Equation [4] (assuming 
n ~ n e ) and setting it equal to Equation [5l we can get a 
relation between the mass-loading parameter K, the Ha 
luminosity, and the radius (Ofek et al. 2013b) 



4:Trhv H a cn (3K Z 
«2xl0 39 M 2 01 «w,600i8»-r5 1 ergs- 1 J 



(6) 



The reason for the inequality is that it is possible that 
not all of the hydrogen is ionized or that the temperature 
of the gas is higher than 10 4 K (i.e., a^ a depends on 
temperature), and that f3 > 1. We note that if j3 < 1 
then this inequality is incorrect. However, it is reasonable 
to assume that the width (i.e., r% — r) of the hydrogen 
envelope is of the same order of magnitude of r. 

Another important caveat (that can be expressed in 
terms of /?) is that unlike the X-ray and radio emission 
which originate at the shock region (radius r), the narrow 
component of the Ha may originate at radii t"h q > r (i.e., 
above the shock region). In this case Lhq cx M^ a r^i, 
then Mn a /M « (rna/r) 1 / 2 . Therefore, if rn a is an order 
of magnitude, or more, larger than r, the contribution of 
Mn a to the bound-free column density (Eq. [2]) and the 
diffusion time scale discussed in £13.41 will be small. 



3.3. X-ray emission 

The X-ray emission from optically thin region is given 
by (e.g., Immlcr et al. 2008) 



4irr 2 A(T)n 2 dr, 



(7) 



where A(T) is the effective cooling function in the 
0.2-10 keV range. Assuming an optically thin thermal 
plasma with a temperature in the range 10 6 -10 8 K (Ray- 
mond ct al. 1976), we adopt a value of A(T) w 3 x 
10 -23 erg cm 3 s _1 . Substituting Equation Q] into Equa- 
tion |7] we get (e.g., Ofek et al. 2013b) 



L x «4ttA(T) 



K 2 



(/z p ) 2 m 2 r 



-(r+Tbf) 



= 3.8 xl0 41 M 2 01 «- 2 00 r^e 



-(r+r bf ) ergs -l^ 



(8) 



This expression includes additional exponential term due 
to absorption in the wind, where r is the Thomson opti- 
cal depth (see Ofek et al. 2010; Balberg & Loeb 2011), 
which is ~ 0.3Mo.oi w w 500 r i5 • Although the Thomson 
optical depth is well known, when the optical depth is 
of the order of a few, Compton scattering is expected 
to reprocess more energetic photons into the 0.2-10 keV 
band (Chevalier & Irwin 2012; Svirski et al. 2012). 
Since the exact X-ray spectrum is not known (Katz et 
al. 2011; Svirski et al. 2012), a proper calculation of 
Lx when r <; 1 is not straightforward. Assuming that 
the Comptonization of hard X-rays into the soft X-rays 
band is smaller than the reduction of soft X-ray flux by 



the optical-depth factors, Equation [5] provides an order 
of magnitude lower limit on M. 

In Figure [U the black-solid line shows an order of mag- 
nitude estimate of the expected X-ray luminosity, assum- 
ing an optically-thin wind-profile CSM with a mass-loss 
rate of 7xlO _4 M Q yr _1 and v w = 500 km s -1 (eqs.[8]and 
fT2]) . Some points are discrepant by factors of 2-3 in lu- 
minosity from this estimate. However, this is a simplistic 

1/2 

model and since M oc LJ , our mass-loss estimate based 
on the X-ray luminosity is plausibly correct to within an 
order of magnitude. 

We conclude that this formula can be trusted only for 
M <y 10 -2 M yr -1 . Above this mass- loss rate, r and Tbf 
are larger than unity. 

3.4. Diffusion time scale 

Another observable that can be used to constrain the 
mass-loss rate is the rise time of the SN light curve. If a 
considerable amount of material is present between the 
SN and the observer, then photon diffusion will slow 
down the rise time of the SN light curve. Therefore, 
the maximum observed SN rise time can be used to put 
an upper limit on the amount of mass between the SN 
and the observer. The diffusion time scale in an infinite 
wind profile is given by e.g., Ginzburg & Balberg (2012) 

kK r , / c \ 

tdm~ In — -1 

c \v sh J 

£* 0.13 K o.34Mo.oifw 1 5oo[ h i(30f s ~h! 4 ) - 1] day. (9) 

Here v s h.4 is the SN shock velocity in units of 10 4 kms^ 1 . 
In the case of SN 2009ip the early fast rise of the SN light 
curve provides an upper limit on tdis and, therefore, an 
upper limit on M. 

3.5. Free-free absorption 

Typically, SN progenitors with mass-loss rates of ~ 
10 -6 MQyr -1 are easily detectable in radio frequencies 
in the nearby universe (e.g., Horesh et al. 2012; Krauss et 
al. 2012). The radio emission is the result of an interac- 
tion between the SN shock and the CSM that generates 
synchrotron radiation peaking at radio frequencies (e.g., 
Slysh 1990; Chevalier & Fransson 1994; Chevalier 1998). 
However, if the material is ionized or partially ionized 
then the free-free optical depth may block this radiation. 
The free-free optical depth in a wind profile, between 
radius r and the observer, is given by (e.g., Ofek et al. 
2013a) 



-3 



1.0 x l^T-l ^ \-% Ml m r t i 



(10) 



where v\q is the frequency in units of 10 GHz. We note 
that the presence of Balmer lines in the spectrum likely 
means that at least some of the Hydrogen is ionized and 
therefore, free-free absorption is important. 

Chandra & Soderberg (2012) and Hancock et al. 
(2012) reported on radio observations of SN2009ip ob- 
tained on 2012 September 26 using the Jansky Very 
Large Array ( JVLAu) , and the Australia Telescope Com- 

9 The Jansky Very Large Array is operated by the National 
Radio Astronomy Observatory (NRAO), a facility of the National 
Science Foundation operated under cooperative agreement by As- 
sociated Universities, Inc. 
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pact Array (ATCA). The JVLA observations did not de- 
tect the SN in the 22 GHz and 8.9 GHz bands down to 
a 3-fi upper limits of 131 /iJy and 65 /iJy, respectively. 
The ATCA observations, put a 3-cr limit of 66 /xJy in the 
18 GHz band. Since our previous limits show that there 
is a significant amount of CSM interacting with the SN 
shock in this event, it is likely that strong synchrotron 
radiation is generated. The non-detection of such a radio 
source implies that rg > 1, providing, therefore, a lower 
limit on M. 



3.6. Precursor fluence 

Prior to the fast rise detected on 2012 September 24 
(Prieto et al. 2012), the light curve of SN2009ip pre- 
sented a feature which can be interpreted as ^ 1 month- 
long outburst. If we assume that this outburst was a 
mass-loss event (rather than part of the SN explosion), 
and if we assume that the bolometric luminosity of the 
outburst is of the same order of magnitude as the kinetic 
energy released in the outburst, then by comparing the 
bolometric luminosity with the kinetic energy we can get 
a rough estimate of the mass released in the outburst. 
The outburst had a peak absolute U-band magnitude of 
about —15, and a duration of at least 30 days (see Pri- 
eto et al. 2012). Therefore, the total bolometric fluence 
of the outburst is Ebol ii 8 x 10 47 erg. The reason for 
the lower limit is that we do not know the outburst light 
curve bolometric correction^, and we only have a lower 
limit on its duration. Comparing -Eboi with the kinetic 
energy and dividing by the duration of the event, idur, 
we get a lower limit on the mass-loss rate 



V Edur 

-1.6 X lO-\- 2 2ooo ^ bo l,8 e 47^ u 1 r,3O M 0y^ 1 J ( U ) 

where u Wj 2ooo is the wind/outburst velocity in units of 
2000 km s -1 , ^01,8647 is the bolometric energy in units 
of 8 x 10 47 erg, and tdur,30 is the outburst duration in 
units of 30 days. 

4. DISCUSSION 

4.1. Constraints on Mass loss 

Along with the observations, equations [2 [6l [8j E3 [10] 
and [TT] provide order of magnitude lower and upper 
bounds on the mass-loss rate from the SN progenitor. 
It is important to note that the reason these are only 
order of magnitude estimates is because some of the as- 
sumptions that go into these formulae are likely inaccu- 
rate. For example, the assumption that the wind is infi- 
nite, continuous and can be described by a single veloc- 
ity component, or the assumption of spherical symmetry. 
Nevertheless, these relations provide order of magnitude, 
independent, estimators for the SN progenitor mass-loss 
rate. 

Some of the estimators require knowledge regarding 
the shock radius r. Following Chevalier (1982), we use 

10 The bolometric magnitude correction is always positive. 
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Fig. 4. — Upper and lower limits on the mass-loss rate of the 
SN 2009ip progenitor as a function of the wind /outburst velocity. 
The assumptions that go into these bounds are discussed in |J3] 
The solid black line shows the limit based on the X-ray luminosity 
(Equation [Sj . The solid gray line represents the Ha luminosity 
derived limits (Equation [B| . The black dashed line and the gray 
dashed line represent the column density (Equation[2]| and diffusion 
time scale (Equation [9]l limits, respectively. The verti cal t hick gray 
solid line is based on the bolometric fluence (Equation l 111) , Finally, 
the free-free absorption limit (Equation HOjl is represented by the 
black dash-dotted line. The arrow attached to each line marks the 
direction of the region allowed by the line criteria. The horizontal 
dotted lines mark the 200 and 2000 km s — 1 wind/outburst velocity 
which we derive from the, presumably pre explosion, IR spectrum 
(Figure [SJ. These lines constitute the approximate range of plau- 
sible wind/outburst velocities. Since it is likely that the mass-loss 
was not a continuous process with constant mass-loss rate, these 
measurements should be regarded as an order of magnitude esti- 
mate. See discussion in [|4] 

the approximation^ 

r ~f v (t)dt = (5/4)u cj , 8 (i 8 - t ) 1/5 (t - io) 4/5 

/ \ 4 /5 

~ 3 x 10 15 "i 3 ' 3 -i ( ^ttt 2 - ) cm, (12) 

8000 km s 1 ^ 30 day J ' v ' 

where t is the time, to is the SN explosion time, and v e j jS 
is the SN ejecta velocity (« 8000 km s^ 1 ) at t s (t s — t = 
71 day), In Figure0]we present the limits we derive on M, 
as a function of the wind/outburst velocity. Specifically, 
assuming N ~ Ah, our X-ray observations of SN2009ip 
provide an upper limit of -ZVh < 2 x 10 22 cm -2 . There- 
fore, Equtaion [5] constitutes an upper limit on M (black 
dashed line). Using Equation [6[ the Ha line luminosity 
of the narrow Ha component we measured on 2012 Dec 
4, Lhq ~ 1-6 x 10 39 ergs _1 , and assuming (3 = 1, we can 
set a lower limit on M (gray solid line). As discussed in 
ij3l Equation [51 along with our measured X-ray luminos- 
ity of 1.1 x 10 39 ergs _1 sets a lower limit on M which is 
shown as the black solid line in Figure [U Furthermore, 
the SN rise time of 2.3 mag day -1 (Prieto et al. 2012) im- 
plies idiff 55 0.3 day. Along with Equation [9] these provide 
an upper limit on M (gray dashed line). The estimate 
based on the bolometric fluence (Equation [TTj) is shown 

11 Assuming the power-law index describing the ejecta velocity 
distribution is m = 8. Note that Chevalier (1982) denoted this 
variable by n, while Balberg & Loeb (2010) use m. 
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as a vertical gray thick line. The non detection in radio 
(Chandra & Soderberg 2012; Hancock et al. 2012) with 

Equation [TO] set a lower limit on M (black dash-dotted 
line). On each line we mark also the (t — to) used to 
calculate the line position. Here we assume that to is on 
2012 Sep 24. 

Although there is no single M value which is consis- 
tent with all the bounds in Figure |4j the closest val- 
ues to all the bounds is in the range of about 10~ 3 to 
10- 2 M Q yr- 1 . 

Kiewe et al. (2012) review the observed properties of 
15 type-IIn SNe. They reported mass-loss rates, prior 
to explosion, in the range of 10~ 4 to ~ IMoyr" 1 , 
while their wind velocities are in the range of ~ 30 
to 1600 km s -1 . The mass- loss rate and wind veloc- 
ity of SN2009ip is consistent with these values. An- 
other SN which shows some similarities with SN 2009ip 
is SN2010mc (PTF lOtel; Ofek et al. 2013). This SN 
showed an outburst about one month prior to its explo- 
sion. We note that the high state of SN 2009ip just prior 
to its fast rise (Prieto et al. 2012) can be interpreted as 
a similar outburst. 

Type-IIn SNe are likely a non-homogeneous class of 
objects arising from multiple mechanisms. It is not 
clear what the best combination of parameters that will 
help us to relate a given type-IIn to specific mecha- 
nism/progenitor are (e.g., luminous blue variable erup- 
tions). However, better mass-loss rates and wind velocity 
measurements for larger samples of type-IIn SN progeni- 
tors, as well as additional cases of pre-explosion outbursts 
can provide the missing link. 

4.2. Interpretation 

With the exception of the mass estimate based on the 
Ha luminosity, the mass-loss estimators in Figure U are 
consistent to an order of magnitude. Specifically, the 
mass-loss lower limit based on the Ha luminosity, and 
assuming ft = 1, is over an order of magnitude above the 
upper limit which is derived from the bound-free absorp- 
tion column-density limit. We note that in the case of 
another source in which a similar analysis was applied 
(SN2010mc/PTF10tel; Ofek 2012; Ofek et al. 2013b), 
the various estimators were consistent. 

A possible explanation for the discrepancy here is that 
some of our basic assumptions are incorrect. Among 
these assumptions, are uniformity of the CSM, spheri- 
cal symmetry, r~ 2 density profile, solar metallicity, and 
ionized (but not fully ionized) CSM. Alternatively, as 
we discussed in ^3.21 it is possible that the Ha emit- 
ting region is further out, above the shock region (e.g., 
|3 > 1). Therefore, it is possible that the Ha line lumi- 
nosity probes a completely different regions of the CSM, 
than the other methods discussed in $3l 

Interestingly, we note that in Figure@J solution allowed 
by all the mass-loss rate estimators which depend on the 
integral of density along the line of sight (marked by 
non-solid lines in Fig. EJ, infer low mass-loss rates, while 
estimators which measure the total emission from an op- 
tically thin volume (marked by solid lines in Fig. [4]), give 
high values for the mass-loss rate. This behavior hints 
that another possible explanation for the discrepancy be- 
tween the various lines in Fig. [4] is that the CSM has an 
aspherical geometry. 



There are two simple geometries that are roughly con- 
sistent with these results. The first simple explanation is 
that the CSM around the SN has a disk geometry, and we 
observe the system from above or below the disk. In this 
case there will be a relatively small amount of interven- 
ing material between the observer and the source, hence 
small value of Nh, and short ijiff • Moreover, in this case 
the total emission (Lx and Lh«) will be larger relative 
to the expectation based on the spherical geometry as- 
sumption and on the values of Nn and idiff- The second 
simple model is that the CSM has a bipolar hourglass- 
like structure. In this case we observe the system from 
the equatorial plane. We note that there are likely other 
possible geometrical solutions, which are more compli- 
cated. 

We conclude that the best explanation for the discrep- 
ancy between the mass-loss estimators is that the Ha 
emission region is above the shock region (or effectively 
/3 ^> 1), or alternatively that the CSM is aspherical. Un- 
fortunately our order of magnitude analysis does not pro- 
vide a way to distinguish between the two scenarios. 

If the Ha emission region is indeed located further 
out, relative to the shock, then an immediate conclu- 
sion is that an order of magnitude estimate to the mass- 
loss rate during the SN precursor is in the range of 
- 10~ 3 M Q yr" 1 to - 10~ 2 M yr" 1 . 

In order to convert the mass-loss rate to an estimate 
of the total mass in the CSM we need to integrate Equa- 
tion [5] out to a specific radius. Here we choose to inte- 
grate the total mass out to radius of 6 x 10 15 cm. The 
reason for this choice is that the relatively abrupt disap- 
pearance of the X-ray flux ~ 70 days after the explosion 
may indicate that the CSM density is falling (faster than 
a wind profile) at a distance of 6 x 10 15 cm (see Equa- 
tion [12]) . Using Equation [5] we find that the total CSM 
mass out to this radius is 

Mcsm ~ 4 x 10- 2 Mo.oi« w ,5oo^i;i5 M . (13) 

Here r^Ei5 is the radius in units of 6 x 10 15 cm. 

4.3. Implications 

Pastorello et al. (2012) and Mauerhan et al. (2012) 
suggested that the outbursts of SN 2009ip are due to pul- 
sational pair instability. However, Woosley et al. (2007) 
predict that the mass-loss in pair instability mass ejec- 
tions would be at least a few solar masses. Unless the 
geometry is highly aspherical, this theoretically predicted 
mass-loss is high relative to our estimate of the total mass 
in the CSM (i.e., ~ 0.1 M Q ). Our mass- loss estimate is of 
the same order of magnitude as the one derived in Ofek et 
al. (2013b) in context of the Quataert & Shiode (2012) 
mechanism. The estimators presented in Figure |4] are 
also in rough agreement with the shell mass of ~ 0.15 M© 
suggested by Soker & Kashi (2012), in context of their 
binary-star merger scenario. 

Interestingly, both Levesque et al. (2012) and Soker & 
Kashi (2012) suggested aspherical models for SN2009ip 
(see also Mauerhan et al. 2012). Levesque et al. (2012) 
argued for a thin disk geometry, while Soker & Kashi 
(2012) suggested a bipolar hourglass-like geometry. How- 
ever, Soker & Kashi (2012) suggested that we are observ- 
ing the system along the polar direction. We note that if 
the discrepancy in Figure [4] is due to asymmetry in the 
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CSM, rather than the radius at which the Ha line is gen- 
erated, than the Soker & Kashi (2012) geometry is not 
consistent with our observations. However, we do not 
claim that our suggested geometries are the only possi- 
ble solutions. Finally, we note that different mass-loss 
estimators have different functional dependencies on r. 
Therefore, additional observations (e.g., radio) can con- 
strain the density profile of the CSM. 
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